Quantum spin Hall effect in 2D topological insulators 

E. B. Sonin 
Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 91904, Israel 

ABSTRACT 

The original motivation of great interest to topological insulators was the hope to observe the quantum spin Hall 
effect. Therefore if a material is in the topological insulator state they frequently call it the quantum spin Hall 
state. However, despite impressive experimental results confirming the existence of the quantum spin Hall state, 
the quantum spin Hall effect has not yet been detected. After a short overview of what was originally suggested 
as the quantum spin Hall effect (quantum spin conductance determined by the topological Chern number) the 

for\' paper analyzes the crucial role of the boundary condition on the observation of the effect and finally discusses 

^3 . whether and how the quantum spin Hall effect could be observed. 

Keywords: Topological insulator, quantum spin Hall effect, quantum charge conductance, quantum spin con- 
ductance, spin accumulation 
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1. INTRODUCTION 



Topological insulators have recently attracted great attention and their study is developing into a new exciting 
area of condensed matter physicsPKED The key feature of topological insulators is the presence of helical edge 
J^ ' states in 2D systems (or surface states in 3D systems). The edge states cross the whole forbidden gap separating 

the valence-band and the conduction-band bulk continua. Helicity of edge states means that electrons with the 
same spin can move only in one direction, which is opposite for two spin directions. As a result, the edge states 
are robust against elastic backscattering, which conserves spin, and the electron transport along edges becomes 
ballistic. The edge state of the 2D topological insulator were experimentally detected in the HgTe quantum welP 
by studying charge transport. It was demonstrated that at the quantum well thickness exceeding the critical 
value 6.3 nm there was an interval of gate voltages where the conductance reaches the quantum conductance 
value independently of the sample width W (see Fig. [T]). This is a clear evidence of the ballistic transport 
through edge states while the main bulk is not conducting. The topological insulators states were also detected 
in 3D materials PES] Though nowadays topological phases are mostly investigated in solids they were already 
well known in superfluid 3 He from 80sPQ 

Originally 2D topological insulators were introduced as systems, in which the quantum spin Hall effect (QSH) 
was expectedJIEl anc | the state of the topological insulator is frequently called the QSH state. However, it is 
necessary to clearly specify the definition of the spin Hall effect in general and the quantum spin Hall effect 
ff\ , in particular. Historically the phenomenon, which is called nowadays spin Hall (SH) effect, was defined as 

£T) • an edge spin accumulation resulting from a bulk spin current transverse to an external electric fieldP^ At 

the sample boundary, which fully reflects all electrons, the spin current induced by the electric field must 
be compensated by a spin diffusion current. Spin diffusion is a dissipative process, which is related with non- 
equilibrium spin accumulation at the edge (or spin orientation in the terminology of Dyakonov and PereP^) . Since 
spin accumulation is more accessible for observation than spin current, usually observation of spin accumulation is 
considered as an experimental confirmation of the SH effect E2E2 But later it was revealed that spin accumulation 
is not an ideal probe of bulk spin currents since spin accumulation is possible without bulk spin current, and, 
vice versa, spin current not necessarily results in spin accumulation (see the review^ and references therein) . 
This is a consequence of non-conserved total spin. In particular, at the equilibrium bulk spin currents cannot 
lead to spin accumulation, since the latter is accompanied by dissipation. Instead spin is transformed to an 
orbital moment at the sample edge (edge torque^). The same is true for bulk spin currents in insulators: since 
an external electric field does not pump energy into the sample (no Joule heating) dissipative spin accumulation 
cannot be supported. 
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Figure 1. Voltage biased topological insulator. (a)Edge states. Wide arrows show spin direction (spin quantization axis 
is not in the plane as in the figure). At the upper edge rightmovers have spin s up while leftmovers have spins down. At 
the lower edge directions of motion are opposite, (b) Voltage variation between electrodes. The voltage varies only near 
electrodes whereas inside the sample there is no electric field because of ballistic transport along edge channels. 



The adjective "quantum" was added to the SH effect in the topological insulators in order to stress that the 
spin conductance (ratio of the bulk transverse spin current to the electric field) was universally determined by 
the topological Chern number for the 2D Brillouin zoneP Such a QSH effect has not yet been detected, since 
the author is not aware of any measurement, which probes anything connected with the bulk spin conductance 
and no experiment up to now is able to detect the predicted quantum spin conductance in the topological phase. 
Nevertheless, it is believed that the QSH effect has already been observed. 4 The ingenious experiments of Konig 
et alW definitely confirmed existence of ballistic edge states and revealed the quantum charge conductance but 
there were no direct or indirect evidences of quantum spin conductance simply because nothing dependent on 
spin was measured. Though both, the ballistic edge states (and the quantum charge conductance as a result 
of it) and the quantum spin conductance, originated from topology, the common origin does not make two 
phenomena identical. An essential progress in studying spin-related phenomena in topological insulators was 
recent detection of spin polarization of the edge states^ using the inverse SH effect, but this new impressive 
achievement still does not provide any information on bulk quantum spin conductance, since the latter does not 
affect any measurement result. So the QSH effect, at least in its original meaning, still waits its experimental 
confirmation. 

The present paper elaborates this point of view. It starts from a short overview of the model by Bernevig, 
Hughes, and Zhang (BHZ) for the topological insulator (Sec. ^ and the definition of spin and spin current in 
this model (Sec. [3]). In Sec. 0] the quantum spin conductance for the periodic boundary is discussed, which 
is determined by the topological Chern number and was originally considered as the determining feature of the 
QSH effect. Section [5] shows that in the strip geometry with fully reflective lateral edges it is impossible to 
realize the quantum spin conductance. But in this geometry there is non-quantum spin conductance, which is 
possible both in the topological and conventional-insulator phases of the BHZ model. Section O addresses spin 
accumulation (polarization) at edge states and argues that observation of spin polarization cannot be considered 
as demonstration of the QSH effect. The concluding Sec. [7] discusses a possibility to detect the genuine QSH effect 
using the Corbino geometry and summarizes the present analysis: the QSH effect is still waiting its experimental 
confirmation. 



2. THE TOPOLOGICAL-INSULATOR MODEL OF BERNEVIG, HUGHES, AND 

ZHANG 

The simplest and therefore the most transparent model for studying fundamental features of the topological 
insulator is the model suggested by Bernevig, Hughes, and Zhang (BHZ) for the topological insulator in the 
HgTe quantum wellPQES The model is a simplified version of the Kane model and its 4x4 Hamiltonian is given 

where 

H(k) = e(k)I + dm (2) 

is a 2x2 Hamiltonian, fi are Pauli matrices of the pseudospin, % = x,y, z, and 7 is a unit 2x2 matrix. Assuming 
that all essential processes occur at low k close to the Brillouin zone center, the components di of the vector d 
are 

d x = Ak x , dy = Ak y , d z = e (fc) = M - Bk 2 . (3) 

Two components of the pseudospin in any 2x2 block of the Hamiltonian ([1]) correspond to the valence (pseudospin 
up) and the conduction (pseudospin down) bands, which overlap in the topological insulator phase at M > 0. 
The off-diagonal linear in k terms in any block lead to mixing of two original bands and to forming new bands 
separated by a forbidden gap. The conventional-insulator phase without edge states corresponds to the condition 
M < 0. Further we shall neglect e(k) in the Hamiltonian @ as not important for the outcome of the analysispJ 

The lower block yields the states obtained from those for the upper block by the time-reversal transformation. 
Because of the absence of off-diagonal blocks in the Hamiltonian ([1]), one can analyze states for any block 
separately. We consider the upper block. The eigenstates in the k space for the Hamiltonian ([2|) are spinors 



1 / Ve±eo \ 
*±(k x ,k v ) = -==[ A (k x+i k v ) , (4) 



where e = \d\ = yeg + A 2 k 2 , and the spinors \&+ and >&_ correspond to the energies +e and — e respectively. 
At M < A 2 1 IB the energy of the upper band has a minimum at k — ( Fig. [2h.). At M > A 2 /2B the 
energy has a maximum at k — 0, whereas the minimum band energy e TO = Ay/M/B — A 2 /AB 2 corresponds to 
k m = y/M/B - A 2 /2B 2 (Fig. [2b). 

An edge state near the edge y = should be a a superposition of two states of the same energy and k x : 



ik vlV _J_ n „\b , (h_ iv n \0 ik yZl>~] P ik * x 



* = [ai*+(fc XI k yl )e ik " iv + a 2 ^ + (k x ,k y2 )e 



e 



(5) 



where k y \ and k y2 are two complex solutions with positive imaginary parts of the bi-quadratic equation for k y 
following from the energy spectrum at fixed e and k x . 

The condition ^ = at the edge y = yields two equations for a\ and a 2 , which have a solution if 

A(k x +ik y i) _ A(k x + ik y2 ) 
e e + e (fci) e e + eo(fe) 

For M > this equation is exactly satisfied if e e = —Ak x , which corresponds to the edge state 

* oc ( J J e**-* (e~ p+y - e~ p - y ) , (7) 



where p± = A/2B ± y/k x — M/B + A 2 /AB 2 . Considering the lower block of the Hamiltonian ((T|) one receives 
the spectrum e e = Ak x . 
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Figure 2. Band energy as a function of k (solid lines) and edge state energy as a function of k x (dotted and dashed lines) 
in a topological insulator (M > 0). Dashed and dotted lines correspond to opposite spin directions [upper and and lower 
block of the Hamiltonian ©]. a) M = 0.3A 2 /B < A 2 /2B; b) M = 2A 2 /B > A 2 /2B. 



3. SPIN AND SPIN CURRENT IN THE BHZ MODEL 

Starting the analysis of spin currents and accumulation, it is necessary to understand what "spin" we would 
like to focus on. In HgTe the conduction band originates^] from a s-type (I = 0) atomic orbital, and its total 
moment coincides with spin, whereas the valence band is related to a p-type (J, — 1) atomic orbital and has the 
total moment j = 3/2 with its projection on the quantization axis (the axis z normal to the insulator plane) 
m,j = 1/2. This projection of the mechanical total angular momentum was also called spin in the previous 
works PEI We shall adopt this extension of the term "spin" remembering that in general this is an effective spin 
but not genuine one. According to this definition the spins in the valence and the conduction bands in HgTe are 
equal to the electron spin Ti/2. However, the mechanical moment would be relevant only if the theory addressed a 
mechanic effect, like that considered in RefP^l If the goal is to describe electromagnetic phenomena like the Kerr 
effect or any other electrodynamic effect, one should use the magnetic moment, which depends on the Lande 
factor of the atomic orbital. The magnetic moments in the the conduction and valence bands of HgTe are fiB 
and 2/xs/3, where [Lb — eh/2mc is the Bohr magneton. In order to take it into account we further generalize 
the definition of effective spin, introducing "spins" s c and s v for the conduction and the valence bands. If the 
analysis addresses a mechanical effect then s c — s v = %/2 in HgTe quantum wells as was supposed earlierPSl 
But if the magnetic moment is relevant for the theory, we define effective spins of the conduction and the valence 
bands as s c = Ti/2 and s v — fi/Z, i.e., as magnetic moments of the bands divided by the electron gyromagnetic 
ratio. So in general the band effective spins s c and s v differ and the operator of the effective spin in the BHS 



model is given by^S 

s z = s z I + As z f z , (8) 

where s z — (s c + s„)/2, and As z = (s v — s c )/2. 

Let us consider now the balance equation for spin (the continuity equation with the torque term) for an 
arbitrary quantum state in the BHS model. The balance equation can be derived from the Schrodinger equation 
as was explained in details for the Rashba HamiltonianP^ Restricting ourselves with the z component of the 
spin density S z = \&^s 2 \I>, the balance equation is 

^± + Va .r a = G\ (9) 

where the torque is 

G z = iAs z A |* f • [V x i-] z * + [V x t] z & ■ *} , (10) 

and the spin current is given by 

J z = -* f {s z ii + Vi§ z } * = s z &Vi<& + As^oi* 1 *. 



Here 

dH(k) 



hdki 
is the group velocity operator and voi(k) — de^/fidki. 

The first term in the spin current Eq. (|lll) differs from the charge current 

J, = e* 1 ^,* 



(11) 
v 0i (k)f z + An (12) 



(13) 



by the constant factor s z /e. Only this term was taken into account in previous publications assuming As 2 = 0. 
In this approximation our effective spin is a conserved quantity and the torque G z vanishes. But in general the 
second term proportional to the difference of spins in two bands should not be ignored and leads to important 
effects as shown below. 

An interesting consequence of helical edges in the topological phase M > is a persistent spin current flowing 
around the sample. According to Eq. (J7J), two edge states transport the average spin ±s z , and the spin current 
along the edge is: 

f = s z (n^ +n^)v e . (14) 

where v e — de e /fidk x = A/Ti is the group velocity at edge states. This current exists even in the equilibriumj221 
when there is no external electric field and the ID densities n_>. and n<- of right-moving and left-moving charge 
carriers are equal. So this is another example of an equilibrium spin current P^ 

4. THE EFFECT OF AN ELECTRIC FIELD: QUANTUM SPIN CONDUCTIVITY 

The QSH effect follows from the calculation of the spin current induced by an electric field using the Kubo ap- 
proachP The approach is based on the perturbation theory. In the presence of the electric field the perturbation- 
theory expression for the quantum-state spinors are 

*±(lC,ky) = e*W« U ± + ^v x * ± X . (15) 
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Figure 3. Mapping of the Brillouin zone on the spherical surface of the unit vector d = d/d. (a) The conventional 
insulator, M < 0. Any path in the Brillouin zone starting in its center and ending at its border maps on the path in the 
d space starting and ending at the same pole of the sphere. The total area subtended by the vector d vanishes, (b) The 
topological insulator, M > 0. Any path in the Brillouin zone starting in its center and ending at its border maps on the 
path in the d space starting at one pole and ending at the other. The vector d subtends the whole spherical surface with 
area 4ir. 



The transverse spin current in this state follows from its definition by Eq. (|11[) and is equal to 



J v \kx, k y ) 



iheE 
4e 2 



& ! ±[v y ,v x ]*± + As z v ay = s z 



eE 

2h~ 



As Z V y: 



where the term 



A 2 
Q = — (e - hkv ) = d ■ 



dd dd 



is responsible for the topological contribution to the current. Its integral 

0,1k = 2a (1 + i^rl = J -<< ' 



(16) 



(17) 



(18) 



is equal to the area of the spherical surface subtended by the unit vector d — d/\d\ when integrating over the 
whole 2D Brillouin zone. Here C is the Chern number, which is equal to 1 in the topological-insulator phase 
M > 0, but vanishes in the conventional-insulator phase M < 0, in full agreement with topological theorems. 
This is illustrated in Fig. [3] 

After integration the single state current (|16[) over the k space and summation of the contributions of the 
two blocks in the Hamiltonian (JXJ) with opposite directions of spin, only the topological term contributes to the 
total spin current J z = a s E. Eventually the spin conductance is 



2es z 



-C. 



(19) 



This derivation of the quantum spin conductance is well knownp but it is valid only for periodical boundary 
conditions. 



5. EFFECT OF REFLECTIVE BORDERS ON THE BULK SPIN CURRENT 

Let us consider now fully reflective edges of the 2D topological insulator. If the spins of two bands are equal, i.e., 
As z — in the expression Eq. (I16p for the spin current, the spin current is proportional to the charge current. 
In the BHZ model spin- flips processes are absent [no off-diagonal blocks in the Hamiltonian Eq. (JTJ], the total 
number of electrons with given spin is conserved, and their current normal to the reflective border should vanish. 
Correspondingly the spin current, which is a counterflow of electrons with opposite spins must also vanish. So 
the bulk spin current is possible only due to the second term ex As z in Eq. (Tl6|) . A proper eigenfunction near 
the reflective border is a superposition of an incident and a reflected wave a\ i &±(k x , k y ) + a,2^±{k x , —k y ). The 
superposition must satisfy the condition \a\\ J y {k x , k y )-\-\a2\ 2 J y {k x , ~k y ) = imposed by the charge conservation 
law. Here J y (k x ,k y ) is the charge current component normal to the border, which is determined by Eq. (|13[1 
for the wave function Eq. (|15|) taking into account the electric field. In the presence of the electric field |ai| 2 
and 1 0,2 1 2 are not equal. According to Eq. (fT6|) and assuming that the electric field does not change the average 
density (|ai| 2 + \a 2 \ 2 = 2), the intensities are 

As a result, the term in the spin current proportional to the average spin s z vanishes but the term proportional 
to the spin difference As z still remains: 

J z (k x ,k y ) = A s 2 « 0y (| ai | 2 - |o, 2 | 2 ) = As z ^g^. (21) 

The total spin current in the whole band does not reduces to the Chern term and is determined by the integral, 
which does not vanish in a conventional insulator (M < 0): LJ 



/ 



A i A 2 -2MB+A^A^AMB , t A 2 -2Mb \ nf u ^ A 2 

G^-dk = 27T I ^ 2VA2-4MB A 2 -2MB-AVA 2 -4MB " t " 2\M\B ) dl Iu *- 4S , , 

J " '' ' / A arctan MMMEM W^ ^-™^ at *L <r M <T £- 

sJTMB~~A? A 2 -2MB +m 2\M\B J dl AB ^ m ^ 2B 



The parameter C for the periodic and reflective boundary conditions is shown in Fig. [¥] Thus, in contrast to the 
analysis based on plane-wave eigenstates satisfying the periodic boundary condition, the bulk spin current may 
appear both in the conventional and the topological insulator, and is not governed by the Chern number if the 
edge of the sample is fully reflective. 

6. BULK SPIN CURRENT AND SPIN ACCUMULATION 

Observation of the calculated bulk spin conductance in the experiment encounters the problem that in an ideal 
topological insulator there is no electric field in the bulk since ballistic edge channels, which are robust against 
elastic scattering, shortcircuit the bulk transforming a topological insulator into a ballistic metal. According to 
the Landauer-Buttiker theory of ballistic conductors the whole voltage bias between electrodes drops only in the 
contacts but not in the bulk (see Fig. [TJb)). Correspondingly the bulk spin current linear in the electric field 
cannot appear and it is impossible to measure quantum spin conductance (i.e., to demonstrate the QSH effect 
in its original definition) in the strip geometry shown in Fig. [1] 

The electric field parallel to the axis x may appear only if there are magnetic impurities at the sample edge, 
which make carrier backscattering possible and lead to finite conductance of edge channels. However in the 
absence of the spin conservation law the bulk current not necessarily leads to accumulation. It may result in an 
edge torque without accumulation, as in the case of equilibrium spin currents in the Rashba medium. 15 As already 
mentioned in Introduction non-equilibrium spin accumulation is impossible also in insulators since the electric 
field does not pump energy into the material without longitudinal charge current. The BHZ model is a model of 



*One cannot extend this analysis on M > A 2 1 IB since in this case there are two real values of k 2 . for the same energy 
(see Fig. [2b). So there are two reflected propagating waves, and the spin current cannot be calculated without a detailed 
analysis of the wave function near the edge. 



conventional insulator 



topological insulator 



C 



2.5 

2.(1 










J 


1.5 


~ : 








/ 


1 fl 










s 


0.5 

on 


- 




^^S 


- 




-0.4 


-0.2 


0.0 0.2 


0.4 



MB/A 2 



Figure 4. The parameter C determining the spin conductance. For the periodic boundary condition C reduces to the 
Chern numbers C — 1 and C — in the topological and the conventional phases respectively (dotted line). For the 
reflective boundary condition C is continuous and does not vanish in the both phases (solid line). 

non-interacting electrons, and the band electrons in the bulk and the electrons in the edge states constitute two 
independent non-interacting ensembles. The bulk band electrons are in the insulating phase and cannot produce 
any accumulation of spin. Thus independently from whether the electric field and the related bulk spin current 
do appear inside the sample or do not, it cannot lead to spin accumulation at the edges. Nevertheless, it does 
not rule out spin accumulation in the edge states, which is inevitably appears due to a rather trivial reason: the 
voltage drop between electrodes generates a ballistic electric current J = et> e (n_> — n<_) in edge channels, which 
leads to inequality of the rightmover and leftmover densities n_> and n<_ . Since rightmovers and leftmovers have 
opposite spins this leads to spin accumulation (or polarization), 



S z 



-J, 



(23) 



totally independent from what is going on with bulk spin currents. Only the average spin appears in this 
expression since according to Eq. ([7]) the edge state is a superposition of two bands with equal amplitudes. 

One can rewrite Eq. (|23[) for the edge spin polarization using the relation connecting the edge channel current 
J = (e 2 /h)V with the voltage between electrodes via the charge quantum conductance e 2 /h: 



s z eV 
h v. 



(24) 



where the half of the quantum spin conductance appears as a factor. One may argue that since spin polarization 
contains the quantum of spin conductance its measurement would be a demonstration of the QSH effect. Such 
conclusion is at least debatable. The conductance quantum appears in this expression because of charge conduc- 
tance of the ballistic transport without any connection with the actual bulk spin current, which is proportional 
to As z whereas the spin polarization is proportional to the average spin s z . The choice of the nomenclature is 



always a matter of taste and convention and one may redefine the QSH effect as the effect of spin polarization. 
But such a definition would be very far from the original concept of the SH effectpl which considered accumu- 
lation of spin brought to the edge by the bulk spin current linear in the bulk electric field. For comparison, let 
us mention that the quantum charge conductance appears not only in the quantum Hall (now we mean not spin 
Hall!) effect, but also in other phenomena, e.g. in the ballistic transport in ID channels (the edge state is an 
example of such a channel). But this does not make these phenomena identical to the quantum Hall effect. 

7. CONCLUSION: OBSERVATION OF THE GENUINE QSH EFFECT 

The analysis presented above leads to the conclusion that the genuine QSH effect (bulk quantum spin conductance 
transverse to an external electric field) cannot be detected not only from transport measurements of charge 
conductance but even from observation of spin polarization at sample edges reported recently!^ Indeed, the 
quantum spin conductance determined by the Chern number is not relevant for this observation and was not 
used for its interpretation. The detection of the QSH effect is possible only if experimentalists are able to 
probe the bulk spin current in the geometry with the periodic boundary conditions when there are no edges 
parallel to the electric field and correspondingly no edge spin polarization . There is no doubt that it is a serious 
challenge for experimentalists. Nevertheless such a geometry can be considered, at least in the theory. Figure [5] 
shows the Corbino geometry for a HgTe quantum well, where the transition to the topological phase is possible. 
Applying a radial electric field one creates an azimuthal spin current. The Corbino geometry provides the 
periodic boundary condition, which allows to observe the quantum spin conductance determined by the Chern 
number. With help of a local gate (see the figure) one can create a potential barrier, which leads to reflection of 
electrons propagating azimuthally and therefore can transform the boundary condition from periodic to purely 
reflective. This would make possible to experimentally check two predictions: quantum spin conductance for 
periodic boundary conditions, i.e., the genuine QSH effect, and non-quantum spin conductance for reflective 
boundary condition, which may exist both in a conventional and a topological insulator. An important feature 
of the Corbino geometry with the periodic boundary condition is that there are no ballistic helical edge channels, 
which are parallel to the electric field and shortcircuit the electrodes. So the electric field is present in the sample. 

The most difficult part is detection of a bulk spin current in the absence of any spin polarization. A possible 
method of bulk moment current detection is observation of an electric field generated by any moving magnetic 
moment E5 Pipz| This is the "inverse spin Hall effect" , which has already been observed in a number of exper- 
imentsPS^Sl But in the geometry of the QSH effect the electric field generated by the spin current is exactly 
parallel to the external electric field. So measurement of the spin current requires detection of a small correction 
to the dielectric permeability. The local gate can help to this detection, since it works as a valve, which regulates 
the spin current and therefore tunes the small correction to the dielectric permeability. 

One may also hope to find suitable optical methods of spin-current detection. Recently a second-order non- 
linear optical effect of spin currents was revealed theoretically 27 and experimentally. 28 It was demonstrated that 
a pure spin current, even without spin polarization, is able to induce a second-harmonic generation. Concluding, 
the experimental detection of the Chern number, which determines the quantum spin conductance, seems elusive 
at the present moment, and some new ingenious set-ups should be looked for this goal. 
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